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I.  mntoDiicTzoH  Ain>  poRmnAnoH 

An  laportant  application  of  tha  wing-body  jtincture  problem  is  in 
the  study  of  the  horseshoe  vortex  which  occurs  in  the  vicinity  of 
a  wing/strut-«rall  intersection  for  Instance.  Nxwerous  flow 
visualization  studies  (Thwaites  1960,  Baker  1979,  Barber  1978, 
Thomas  1987)  have  shown  that  the  horseshoe  vortex  system  consists 
of  a  three-dimensional  boundary-layer  separation  in  front  of  the 
wing/strut  followed  by  a  vortex  flow  which  wraps  around  the 
structure.  This  vortical  flow  is  characterized  by  three- 
dimensionality,  unsteadiness,  and  large  Reynolds  number  but 
having  viscous  effects  which  are  crucial  to  the  flow  evolution. 
Such  flows  exist  in  many  situations.  For  example,  horseshoe 
vortex  flow  occurs  near  the  juiKrtion  of  an  airplane  wing  with  the 
fuselage  (Thwaites  1960)  or  the  junction  of  plate  and  support  in 
a  plate  heat  exchanger.  Another  example  is  present  in  axial 
turbomachinery  (Barber  1978)  where  boundary  layers  idiich  develop 
on  the  annular  surfaces  of  the  axial  flow  passage  encounter  rows 
of  stationary  and  rotating  blades.  The  horseshoe  vortical  flow 
is  of  engineering  interest  because  it  can  lead  to  flow 
degradation,  high  wall  shear  stresses,  and  high  local  heat 
transfer  rates.  It  also  plays  a  role  in  the  origin  of  three- 
dimensional  comer  flows. 

Three-dimensional  fluid  motion  past  a  wing-body  juncture  or 
the  related  problem  of  comer  flow  is  of  fundamental  importance 
in  fluid  mechanics  as  noted  above.  In  all  of  these  cases,  the 
concern  is  with  both  the  local  and  global  scales  of  the  viscous 


forcM  produced  by  the  body  as  tiell  as  with  any  significant 
secondary  flows  which  say  be  set  up.  It  is  hoped  that  an 
understanding  of  these  types  of  flows  and  their  scales  will 
provide  us  with  a  sore  cosplete  understanding  of  cosplicated 
three>^isensional ,  high  Reynolds  nuaber  flows.  To  this  end,  the 
concepts  of  interaction  theory  have  been  utilized  with  the 
ultiaate  goal  of  elucidating  the  aajor  effects  of  viscosity  and 
heat  transfer  on  the  local  flowfield  near  a  wing-body  juncture 
for  subsonic  aainstreaa  flows. 

Saith  &  Gajjar  (1984)  were  the  first  to  consider  the  problem 
of  a  steady,  incoapressible,  laainar  flow  past  a  juncture  aade  up 
of  a  localized  saall-scale  thin  wing  protruding  froa  a  locally 
flat  surface  when  the  Reynolds  nuaber.  Re,  is  large,  using  the 
triple-deck  theory  of  viscous-inviscid  Interactions.  A  schematic 
of  the  geometry  considered  by  these  authors  and  that  to  be 
considered  in  this  report  is  shmm  in  Figure  1.  Saith  &  Gaj jar's 
theory  was  developed  to  accomaodate  small  wings  of  finite  span 
which  scale  with  the  streaawise  and  spanwise  triple-deck  length 
scales.  However,  as  was  deterained  as  part  of  this  current 
research  effort,  Saith  &  6a j jar's  theory  is  not  generally 
applicable  to  all  suitably  scaled  wing  shapes,  as  will  be 
discussed  below.  The  incident  boundary  layer  on  the  body  surface 
is  assumed  to  be  driven  by  a  locally  uniform  external  flow  and  it 
is  taken  to  be  well  developed,  attached,  «uid  planar  in  character, 
under  these  conditions  an  interaction  between  the  body  siarface 
flow  and  the  wing  is  established  tihich  presumably  has  the  triple- 
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deck  structuTtt  wherein  the  eotlon  past  the  thin  wing  provokes  on 
the  body  surface  an  unknown  pressure  force  which  interacts  with 
the  unknown  displaceaent  thickness  there.  Due  to  this 
interaction,  the  boundary  layer  on  the  body  surface  becoees 
three^diaensional  and  exhibits  upstrean  influence.  Furtherwore, 
a  nonlinear  Interaction  is  possible  if  the  slope  of  the  thin  wing 
in  the  interaction  region  is  of  the  order  Re'^^*. 

As  shown  by  Smith  &  6a j jar  (1984)  the  fundamental  problem, 
as  sketched  in  Figure  1,  can  be  reduced  to  a  consideration  of  the 
following  boundary-value  problem  for  an  incompressible  steady 
flow  in  terms  of  suitably  scaled  velocity  components  (U,V,W), 
pressure,  P(X,y),  and  displacement  thickness  function,  f(X,Y),  in 
the  (X,Y,Z)  directions,  respectively.  More  complete  details  of 
the  scaling  laws  and  the  triple-deck  nature  of  the  flowfield 
structure  can  be  found  in  Lee  (1994).  The  fundamental  governing 
equations  appropriate  for  the  viscous  layer  in  the  usual  triple- 
deck  theory  are  given  by 


ax  dr  az 


(1) 


ax  dY  dz  ax  az* 


ax  ar  az  ar  az* 


In  addition,  the  following  boundary  conditions  apply 


(3) 


u«v«w-o  onZ«0, 


(4) 
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0  -*  Z  -  <((X,y)t  V  -  0  a*  Z  -»  •,  (5) 

(U,V,lf,P)  -*  (Y, 0,0,0)  as  |X|  -  -,  (6) 

and,  and  according  to  Saith  &  Gajjar  (1984),  the  syaoMtry 
condition 

V  -  0  on  Y  -  0.  (7) 

Tha  problaa  is  closad  sathasatically  by  applying  the 
pressure-displaceaent  thickness  interaction  law,  viz. 


(x-c)*+r* 


* 


_Lr  r  a  Ml 


(8) 


for  two-diaonslonal  wings,  or 


nxt)  -  - 

q^2  ^(X-5)»  +  ira+e* 

1  r*  r* 


(9) 


for  three-diaensional  wings,  where  hf(X)  is  the  scaled  local  two- 
diaensional  wing  surface  shape,  imile  the  three-diaensional  wing 
surface  is  defined  by  hg(X,Y).  Ife  should  note  that  it  is  the 
wing  shape  which  generates  the  forcing  trtiich  in  turn  generates 
the  interaction  and  three-diaensionality  in  the  body's  boundary 
layer.  This  forcing  coaes  froa  outside  the  boundary  layer  itself 
and  is,  therefore,  quite  different  froa  the  usual  three- 
diaensional  interactive  flow  produced  by  a  saall  hui^>  in  a 
boundary-layer  flow. 
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In  general  the  problen  posed  in  equations  (l)  -  (7)  euid  (8) 
or  (9)  above  is  nonlinear  and  thus  requires  a  nueerical 
treatment.  However,  Smith  &  Gajjar  (1984)  and  Lee  (1994)  have 
considered  linearized  solutions  of  these  equations,  using  Fourier 
transform  methods,  valid  for  small  values  of  the  wing  thickness 
parameter  h,  and  several  choices  for  the  wing  shapes  f(X)  and 
g(X,Y). 

These  linearized  results  always  showed  two  maximum  local 
deficits  occurring  in  the  streamwise  surface  shear  stress  on  the 
body,  just  ahead  of  and  beyond  the  front  of  the  wing,  on  the  axis 
of  symmetry,  tdiich  suggests  that  a  regular  separation  or  flow 
reversal  tfould  tend  to  occur  there  first,  when  the  flow  is 
sufficiently  nonlinear.  This  result,  coupled  with  the 
corresponding  spanwise  svurface  shear  stress  patterns,  led  to 
qualitative  agreement  with  experimental  observations,  as  noted  by 
Smith  &  Gajjar  (1984),  of  a  separation  line  starting  upstream  of 
the  wing  and  bending  around  it,  and  of  separation  downstream. 
Furthermore,  the  secondary  flow  patterns  found  from  the 
linearized  analysis  tend  to  tie  in  with  those  found  in 
experiments  also. 

In  view  of  the  encouraging  results  of  the  linearized  theory, 
in  this  research  effort  numerical  solutions  of  the  nonlinear 
boundary~value  problem  defined  by  equations  (1)  ->  (7),  and  (8)  or 
(9)  have  been  found  for  several  wing  shapes  with  the  goal  of 
determining,  among  other  things,  the  conditions  for  boundary- 
layer  separation.  He  have  also  extended  the  theory  to  Include  the 
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•ffects  of  toqporaturo  on  tho  flow  i«h«n  the  nainstrean  is 
subsonic  in  character.  Referring  to  Lee  (1994)  for  details,  it 
can  be  shown  that  the  energy  equation  in  suitably  scaled 
variables  for  the  nondisensional  tesperature  distribution*,  T,  in 
the  viscous  sublayer  is  given  by 


-  JLiL 

^3x  ^ar  ^az  Praz^ 


(10) 


along  with  the  boundary  conditions 

nX,YJ^-Z  as  IKI.IKI  -- 
r(jra^,0)  - o,  nx;r;z)  - z  -  forz-^ 

when  we  can  asstme  that  the  wall  teaperature  and  Prandtl  nuaber, 
Pr,  are  constants.  It  should  also  be  noted  that  the  above  system 
of  equations  is  valid  for  any  svtbsonic  flow  since  all  paraaeters 

related  to  the  coapressibility  effects  of  the  flow  can  be  scaled 
out  in  the  triple-deck  foraulation  of  the  problea  (see  Lee  1994). 

In  obtaining  nuaerical  solutions  to  the  linearized  fora  of 
the  above  equations  for  h  «  1,  Saith  &  6a j  jar  (1984)  introduced 
the  concept  of  an  "effective"  huap  shape  to  transfora  the 
linearized  physical  wing-body  juncture  problea  into  an  equivalent 
three-diaensional  huap-like  problea  using  an  appropriate 
transforaation  relating  the  displaceaent  thickness,  interactive 
pressure  and  the  wing  shape.  This  change  allowed  for  rather 
simple  computations  to  be  done  using  a  finite-  Fourier  transform 
aethod  to  convert  the  closed-fora  solutions  for  the  pressure. 


displaceaent  thickness  and  wall  shear  stresses  in  spectral  space 


7 

to  oqui valent  functions  in  physical  space.  Lee  (1994)  confiraed 
these  results  and  also  coaputed  the  wall  heat  transfer 
distribution  functions. 

In  our  study  of  the  linearized  version  of  this  problea  it 
was  deterained  that  the  theory  developed  by  Saith  &  Gajjar  is  not 
generally  applicable  for  the  airfoils  they  used  in  their  paper 
due  to  the  laposition  of  syaaetry  conditions  on  the  pressure  and 
the  spanwlse  velocity.  These  conditions  are  inconsistent  with  the 
physical  foraulation  of  the  wing-body  problea  except  for  "wing" 
shapes,  f(X),  that  are  coaposed  of  linear  eleaents,  such  as  in  a 
diaaond  shaped  airfoil.  This  restriction  places  rather  serious 
Halts  on  the  types  of  airfoils  for  which  this  theory  is  useful. 
It  should  be  noted  that  atteapts  at  generalizing  this  theory  to 
include  a  broader  class  of  wing  shapes  was  not  successful. 

Of  prlaary  concern  in  this  research  effort  was  the  niiaerlcal 
solution  of  the  nonlinear  governing  equations  (1)  -  (11)  for 
acceptable  wing  shapes  and  order-one  values  of  the  wing  shape 
parameter  h.  Two  different  approaches  for  the  nuaerical  scheae 
were  considered.  In  one  approach,  the  finite-difference  aethod 
for  three-diaenslonal  interacting  flows  developed  by  Bodonyl  & 
DucK  (1988)  irtierein  the  pressure-displaceaent  law,  equation  (8) 
or  ( 9 ) ,  is  replaced  by  a  nuaerical  solution  of  the  upper-deck 
boundary-value  problea  which  is  solved  siaultaneously,  following 
the  ideas  originally  proposed  by  Veldaan  (1979),  along  with  the 
finite-difference  fora  of  equations  (1)  -  (3).  In  previous 
studies  this  approach  has  proved  to  be  quite  useful  in  coaputlng 
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thre«-<lla«n8ional  nonlinear  Interacting  flows.  In  the  second 
approach  the  solution  of  the  system  (1)  -  (3)  along  with  a 
numerical  implementation  of  the  Cauchy-Hilbert  integral,  (8)  or 
(9),  using  Smith's  approach  (1991)  to  discretize  the  pressure- 
displacement,  law  was  determined. 

In  either  approach  the  above  set  of  equations  is  elliptic  in 
the  pressure  field  and  therefore  requires  a  modification  of  the 
multi-sweep  fonrard-marching  methods  which  i«ere  successfully 
developed  for  2-D  problems.  As  noted  by  Smith  (1983)  this 
elliptlcity  causes  an  explosive  3-D  free  interaction  in  any 
forward-marching  procedure  unless  treated  properly.  To  overcome 
these  difficulties,  a  skewed-shear  method  is  used  to  cast  the 
problem  in  a  quasi-two-dimensional  form  (see  Bodonyi  &  Duck  1988 
and  Lee  1994)  while  at  the  same  time  capturing  the  elliptic 
nature  of  the  Inverse  3-D  boundary-layer  problem.  In  terms  of 
these  new  variables,  a  second-order  finite-difference  method  was 
developed  for  the  governing  equations  emd  numerical  solutions  of 
the  finite-difference  equations  were  computed  for  a  variety  of 
grid  sizes  and  several  choic^d  for  the  wing  shape.  The 
computations  are  involved,  requiring  the  use  of  the  CRAY  yilP/864 
supercomputer  at  the  (%io  Supercomputer  Center.  Results  of  the 
numerical  computations  are  given  in  the  following  section. 
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II.  mniSRICAL  RESULTS 

After  considerable  study,  it  was  detersined  that  the  second 
nuaerical  approach  discussed  above  resulted  in  the  aost  accurate 
solutions  of  the  governing  equations.  Therefore,  the  nuaerical 
results  to  be  presented  below  will  be  taXen  froa  that  approach. 
As  noted  above,  for  the  theory  to  be  strictly  applicable,  the 
t%K)-diaensional  wing  shape,  f(X),  aust  be  coaposed  of  linear 
eleaents.  Thus  we  shall  consider  the  diaaond^shapcKl  airfoil 
defined  by 


X  JbrO<X<03 
1-X  ybr  OJ  <  X  <  1 


(12) 


The  grid  sizes  used  in  all  cases  have  lOl  points  for  -4.5  S  X  £ 
5.5,  31  points  for  0  ^  Y  51  3.75  and  41  points  for  0  Z  s  6  in 
the  streaawise,  spanwise  and  noraal  directions,  respectively. 

The  local  convergence  criterion  for  the  interactive  pressure  and 
the  global  convergence  criterion  for  the  velocity  coaponents  are 
10~*  and  10~*,  respectively. 

A  surface  plot  of  the  effective  huiq>  shape  generated  by  the 
diaaond  shaped  airfoil  is  given  in  Figure  2.  It  is  clearly 
evident  that  this  "effective"  huap  shape  is  very  peaked  near  Y«0 
and  X  «  1/2.  This  sharp  peak  places  severe  lialtations  on  idiat 
grid  sizes  are  needed  to  adequately  resolve  the  flow  structure  in 
this  vicinity,  and  it  is  a  prlaary  reason  that  nuaerical 
solutions  could  be  found  only  for  rather  United  ranges  of  the 
paraneters  Involved.  This  situation  is  quite  different  than  in 
the  "regular"  3-D  huap  problen  considered  by  others  since  in 
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tho8«  casM  tha  huap  shape  could  be  aore  conveniently  defined  so 
that  such  situations  as  that  noted  above  did  not  occur.  In  the 
current  study  the  "huap”  shape  is  predeterained,  depending  only 
on  the  actual  wing  shape  and  the  aatheaatlcal  transforaatlon 
inherent  in  the  problea. 

Flgxires  3  through  6  give  coaparisons  between  the  linearized 
solutions  using  the  PFT  aethod  and  the  nonlinear  solutions  using 
the  finite-difference  aethod  for  saall  h.  Note  that  the  results 
are  quite  close  to  each  other,  indicating  that  the  nonlinear 
finite-difference  coaputatlonal  aethod  developed  in  this  effort 
does  %iork:.  For  coapleteness ,  Figures  7-11  give  the  surface 
plots  for  the  pressure,  streaa  and  spanwlse  shear  stress 
distributions,  displaceaent  thickness,  and  wall  heat  transfer 
distributions,  obtained  froa  the  nonlinear  coaputations  for 
h-o.l. 

Nonlinear  solutions  on  the  plane  of  syaaetry,  Y  «  0,  for 
various  values  of  h  are  given  next  in  Figures  12  -  15.  Note  that 
oscillations  begin  to  appear  in  the  solution  for  the  pressure 
when  h  0.5,  and  they  are  clearly  evident  when  h  *  1.0.  This 
nuaerical  instability  could  not  be  reaoved  froa  the  cOTputations, 
and  this  eventually  prevented  the  coaputation  of  further 
solutions  for  larger  values  of  h.  Unfortunately,  this  nuaerical 
breakdown  always  occurred  before  separation  on  the  body  surface 
in  front  of  the  wing  occurred.  Also  note  that  the  nonlinear 
solutions  at  Y  >■  0  do  not  return  to  either  0  or  constant  values 
at  the  downstreaa  liait  of  the  coaputatlonal  doaain.  The  reason 


11 


for  this  discrepancy  is  that  the  cosputational  dosain  for  the 
nonlinear  problea  is  much  saaller  in  extent  than  that  for  the 
linearized  solutions  found  by  the  FFT  eethod.  Due  to  «>eputer 
resource  liaitations  it  was  not  feasible  to  solve  the  nonlinear 
problea  on  larger  doaains,  althouc^  it  is  believed  that  the 
computations  tiould  return  to  their  appropriate  values  if  a  larger 
domain  were  utilized. 

Figure  16  shows  the  flow  streamlines  in  the  plane  of 
symmetry  (Y  «  0)  and  Figures  17  through  19  show  the  cross-flow 
streamlines  near  the  leading  edge,  at  mid-chord,  and  near  the 
trailing  edge  for  h  >  0.5.  Note  that  the  directi  of  the  normal 
velocity  is  downward  near  the  wing  (0  <  X  <  1).  This  effect  is 
indicated  by  the  solid  lines  in  Figures  17  and  18.  However,  the 
normal  velocity  changes  sign  near  the  trailing  edge  as  shown  by 
the  dotted  lines  in  Figure  19.  The  surface  stress  pattern  for 
this  problem  is  shown  in  Figure  20.  The  only  noticeable  effects 
on  the  surface  due  to  the  wing  are  in  the  vicinity  of  the  wing  (0 
<  X  <  1,  Y  -  0). 

The  sharp  peak  of  the  effective  hump  at  Y  •  0  seems  to  be 
the  main  cause  for  the  numerical  difficulties  irtiich  occur  behind 
the  wing.  This  sharp  peak  is  the  result  of  discontinuities  in  the 
derivative  with  respect  to  X  of  the  wing-shape  function  irtiich 
appears  in  the  pressure-displacement  thickness  relationship.  As 
noted  without  reaching  large  enough  values  of  h  it  is  impossible 
to  have  separation  around  the  wing.  For  exaig>le,  the  flow  over  a 
"regular"  hump  shape  h*exp(-x*-y*) ,  will  result  in  boundary-layer 
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Separation  irtian  h  is  approxiaataly  2.8  (saa  Bodonyi  &  Duck  1988). 
Separation  in  front  of  this  wing-body  junction  would  be  axpactad 
to  occur  for  such  order-one  values  of  h.  However,  in  this 
problea,  the  nuaerical  scheae  becoaes  unreliable  for  h  >  0.5,  and 
for  these,  or  saaller,  values  of  h,  the  flow  near  the  surface  is 
still  strongly  attached,  precluding  the  possibility  of  bounr 
layer  separation. 

For  a  study  of  step-size  effects,  an  effective  huap  shape  is 
coaputed  for  a  finer  grid  with  AX  «  0.05  and  AY  «  0.0625  for  the 
dlaaond-shaped  airfoil.  Figure  21  shows  the  coaparison  of  the 
effective  huap  shape  for  the  two  different  aeshes.  The  effective 
huap  shape  generated  using  this  finer  grid  is  input  into  the 
nonlinear  wing-body  solver  using  the  grids  AX  »  0.1  and  AY  • 
0.125.  In  this  case,  the  pressure  shows  oscillations  for  even 
saaller  values  of  h  than  before,  as  shown  in  Figure  22.  Much 
finer  grids  aay  therefore  be  needed  along  with  wider  ranges  for 
each  direction  to  Increase  h  to  larger  values. 

Linear  and  nonlinear  solutions  for  a  three-diaensional  wing 
shapes  of  the  following  foms  have  also  been  coaputed: 

(13) 


8(X.Z.) 


/(J-Z,iai4)#''^  irX>Z,tm4 
0  rX<Z.tSB4 


Here  f(X)  is  the  base  shape  of  the  wing  adjoining  the  surface  and 
Zu  is  the  upper-deck  coordinate,  0  is  a  sweep  angle  and  a  is  a 
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posltiv*  constant  irtiich  controls  tha  dsgrss  of  tapsring  in  the  Z 
direction  avay  from  the  body  junction.  The  results  of  the 
coaqputatlons  are  qualitatively  slsllar  to  those  of  the  two-' 
disenslonal  wings  discussed  above  and  will  not  be  repeated  here. 
For  further  details  on  these  three^dlsensional  solutions  see  Lee 
(1994). 


In  this  r«March  affort  wa  hava  conaldarad  tha  flow 
atructnra  naar  a  snall-scala  wing-body  coabinatlon  within  tha 
fraaawork  of  tripla-daclc  thaory.  Thin  airfoil  thaory  was  usad  to 
obtain  tha  prassura  distribution  around  a  wing  idiich  in  turn 
triggars  a  viscous-inviscid  intaraction  naar  tha  wing-body 
junctura.  As  part  of  tha  formulation  of  tha  problas  wa  hava 
followad  tha  laad  of  smith  &  Gajjar  (1984),  utilizing  tha  concept 
of  an  "affactiva  hui^>  shape"  in  tha  formulation  of  tha  nonlinear 
problem.  This  technique  not  only  simplifies  tha  prassura 
expression  but  also  enhances  tha  convergence  of  tha  numerical 
schema  even  though  tha  concept  itself  is  just  a  transformation  to 
convert  tha  wing-body  problem  into  a  more  conventional  problem 
for  conqnitational  efficiency. 

As  noted  earlier,  there  is  an  important  inconsistency  in 
Smith  8  Gajjar's  paper  on  tha  wing-body  problem.  According  to  the 
triple-deck  scalings,  tha  wing's  presence  enters  in  the  upper 
deck  of  the  triple-deck  theory  as  a  thin-airfoil  formulation  to 
leading  order.  In  the  lower-deck  region,  however,  the  only  effect 
of  the  wing's  presence  is  felt  in  the  pressure  field,  nius  it  can 
be  shown  that  the  boundary  condition  on  the  pressure  must  be 
given  by 

(.» 

However,  Smith  8  Gajjar  omitted  this  condition  and  used  dP/dY  »  0 
instead.  This  is  equivalent  to  taking  the  spanwise  velocity,  V, 
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•qual  to  zoro  at  the  lina  of  syaoMtry,  and  this  is  inconsistent 
for  this  problsB.  As  a  result  it  can  be  proved  froa  the 
linearized  problem  that  there  can  be  no  viscous-inviscid 
interaction  for  general  tiio*^iMnsional  wing  shapes.  Furtheraore, 
this  result  seeas  to  carry  over  to  the  nonlinear  computations. 

The  exception  to  this  result  occurs  for  wings  composed  of  linear 
eleaents  since  in  this  case  f"(X)  is  identically  zero,  lie  note 
here  that  there  can  be  a  viscous-inviscid  interaction  for  a 
three-diaensional  wing  shape  as  discussed  in  Lee  ( 1994 ) . 

It  is  not  clear  at  this  point  idiether  the  triple-deck 
structure  can  accoaaodate  interactions  which  originate  outside 
the  lower-deck  region  as  originally  noted  by  Stewartson  (1974). 
Indeed «  the  results  of  this  study  suggest  that  the  wing-body 
juncture  problem  is  generally  not  of  the  triple-deck  kind 
proposed  by  Smith  &  6a j jar,  due  to  the  inconsistencies  in 
matching  the  solutions  between  the  three  regions  of  the  flow 
field  normal  to  the  surface.  In  particular,  the  effects  of  normal 
pressure  cannot  be  Ignored,  and  the  wing's  presence  must  cone 
into  play  in  a  more  significant  way  than  originally  proposed. 
Unfortunately,  atte^pta  at  developing  such  a  theory  have  not  been 
successful. 
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Figure  1:  Tbe  triple^eck  "box"  describing  the  three-* 

dieensional  flow  about  the  entire  junction  of 
a  thin  sysaetric  wing  is  shown  by  dashed 
lines,  nie  box's  disensions  are  0(Re~*''*)  x 
©(Re-*"*)  X  ©(Re-*"*)  in  the  X,Y,  and  Z 
directions,  respectively,  and  the  distance  to 
the  box  fros  the  leading  of  the  flat  plate  is 
L  -  0(1). 
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Figure  2:  The  effective  hump  shape  generated  by  the  two- 
dimensional  wing  for  a  diamond-shaped  airfoil. 
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Figure  3:  Coaparison  of  surface  pressures  for  linearized  and 
nonlinear  solutions  of  the  t«o-^ieensional  wing- 
body  for  a  diaeond-shaped  airfoil. 
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Figure  4:  Coaparison  of  streamls*  shear  stress  for 
linearized  and  nonlinear  solutions  of  the 
t«fo-diMnsional  wing-body  for  a  diaaond- 
shaped  airfoil. 
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Figur*  6.  COTparison  of  wall  heat  transfar  for  linearized 
and  nonlinear  solutions  of  the  ttfo-diaensional 
wing-body  for  a  diaaond-shaped  airfoil. 
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Figure  7:  Surface  pressure  distribution  for  nonlinear  flow 
(h  =  0.1)  of  the  two-dimensional  wing-body  for  a 
diamond-shaped  airfoil. 
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Figure  8:  Streamwise  shear  stress  distribution  for  nonlinear 
flow  (h  =  0.1)  of  the  two-dimensional  wing-body 
for  a  diamond-shaped  airfoil. 
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Figure  9:  Spanwise  shear  stress  distribution  for  nonlinear 
flow  (h  =  0.1)  of  the  two-dimensional  wing-body 
for  a  diamond-shaped  airfoil. 
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Figure  10:  Displacement  thickness  distribution  for  nonlinear 
flow  (h  =  0.1)  of  the  two-dimensional  wing-body 
for  a  diamond-shaped  airfoil. 
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Figure  11:  Wall  heat  transfer  distribution  for  nonlinear 

flow  (h  =  0.1)  of  the  two-dimensional  wing-body 
for  a  diamond-shaped  airfoil. 
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Figure  12:  Ooapariaon  of  surfaco  prassuras  for  nonlinaar 

solutions  of  tba  two-disansional  wing-body  for  a 
diaaond-sliafMd  airfoil. 


Figure  14:  Coaparison  of  displacownt  thickiMM  for 

nonlinoar  solutions  of  tbs  two-dlasnslonal  wing- 
body  for  a  diaaond-shapsd  airfoil. 


Figure  16:  Flow  strwaalinws  in  thw  piano  of  synetry  (Y  «  0 
for  tho  two-dinonsional  wing-body  for  a  diaaond- 
•hapod  airfoil  (  h  -  0.5,  0  <  X  <  1).  Tho  dottod 
linos  indicato  tho  location  of  tho  wing.  Tho 
contour  valuos  aro  0  (z  «  0),  0.5,  2.0,  4.0,  6.0 
10.0,12.0,  16.0  and  18.0. 


Pigur*  18:  Cross~flow  streaalinas  at:  tbe  mid-cord  (X  « 

0.5)  of  the  two-dimensional  wing  for  a  diamond¬ 
shaped  airfoil  (h  «  0.5).  The  contour  values  are 
-0.0025,0(*),  0.0005,  and  0.01. 


Figure  19:  Cross-flow  strsaalines  near  the  trailing  edge  (X 
•  1)  of  the  two-diaensional  wing  for  a  diamond- 
shaped  airfoil  (h  «  0.5).  The  contour  values  are 
-0.003,  -0.02,  -0.01  and  0(*). 
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Figure  20;  Surface  stress  pattern  for  the  two-dimensional 

wing-body  for  the  diamond-shaped  airfoil  (h=0.5, 
0  <  X  <  1). 


